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SENSORLESS CONTROL OF A HARMONICALLY DRIVEN 
ELECTRODYNAMIC MACHINE FOR A THERMOACOUSTIC DEVICE OR 

VARIABLE LOAD 

REFERENCE TO RELATED APPLICATIONS 
This patent application claims benefit from U.S. provisional patent application 
Serial No. 60/426,048, filed November 12, 2002, the entire contents of which are 
incorporated herein by reference. 

FIELD OF THE INVENTION 
The present invention relates to control techniques for reciprocating linear 
electrodynamic drive machines and, more particularly, to a method of sensorless control 
of a reciprocating linear electrodynamic drive machine for a thermoacoustic device. 

BACKGROUND OF THE INVENTION 

While it is known that sensorless control of rotary electrodynamic machines is 
possible, a number of novel linear machines and relatively new applications for the same 
have motivated the need for robust and cost-effective control methods. 

One such new application where novel and effective control methods for linear 
machines is desired relates to thermoacoustic devices such as thermoacoustic 
refrigerators or similar applications where it may be desirable to track an acoustic 
resonance of the refrigerator for maintaining an optimum operating condition at the 
reciprocating linear electrodynamic machine. 

It is appreciated that many thermoacoustic devices are designed to be co-resonant 
with the linear machine that drives the device. In other words, the thermoacoustic device 
is designed to perform optimally when the mechanical resonant frequency of the linear 
machine is equal to the resonant frequency of the thermoacoustic device. In practice, the 
acoustic resonance can be a function of the operating temperatures or the cooling load, 
and it is then desirable for the excitation provided to the machine, to be adjusted in a 

1 



PST-10902/36 
31011jal 



manner to re-optimize performance. Thus, knowledge of this functional relationship 
between the drive machine and the load device is required for one to control the 
performance of the combined system by controlling the excitation administered to the 
drive machine. 

In conventional thermoacoustic refrigerators, for example, acoustic pressure and 
actuator displacement at the linear machine are typically monitored while a nominally 
harmonic drive frequency and amplitude (either voltage or current) is adjusted to achieve 
a desired operating point in view of these parameters. This desired operating point may 
illustratively correspond to an acoustic resonance at the load, maximum electroacoustic 
conversion efficiency between the linear machine and the load, maximum power delivery 
to the load, or to achieve some other desirable operating point for the thermoacoustic 
system. 

Most often, pressure and displacement are monitored directly using sensors. 
However, U.S. Patent No. 5,342,176 issued to Redlich discloses a method for obtaining 
the piston displacement in a free piston compressor for a linear electrodynamic machine 
using only the measured terminal current and voltage of the machine. Redlich's objective 
for obtaining the piston displacement through the disclosed method was primarily 
concerned with the control of the amplitude of the actuator in the linear machine. Such 
control is necessary to avoid collision with stationary parts of the linear machine, but not 
as a precursory step for controlling a thermoacoustic load, since the application is for use 
with a compressor where efficient compression requires as small as possible a space 
between the piston and the end of the compression space. Further, in the absence of a 
pressure or force determination, it is not possible to determine conditions such as an 
acoustic resonance, the machine efficiency, or to adjust for maximum power transfer. 

U.S. Patent No. 6,289,680 issued to Oh et al. discloses a method of controlling a 
linear motor used in a compressor application wherein force and displacement of the 
linear motor is mentioned, however, there is no suggestion that force or pressure may be 
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extracted from the linear motor in a completely sensorless manner, and these dynamical 
parameters be used to provide a more robust and cost-effective control means. 

Accordingly, there is a need for a method of sensorless control of a linear 
reciprocating electrodynamic machine that obviates the need for pressure and 
displacement transducers traditionally associated with the control of such machines that 
are commonly used for driving thermoacoustic refrigerators or similar frequency 
dependent loads. 

SUMMARY OF THE INVENTION 

The present invention provides a method of sensorless control of a device 
including a linear reciprocating electrodynamic machine used for driving a frequency 
dependent load, wherein sensorless control is accomplished by estimating dynamical 
parameters of the device through the use of a system model and thereafter manipulating 
at least one input to the linear machine to achieve a desired operating condition in view of 
the estimated dynamical parameters and the machine excitation. 

The method for sensorless control begins with providing a thermoacoustic device 
including a linear electrodynamic machine communicating with a load device. The 
inputs to the machine include either the input current or voltage magnitude, and the input 
frequency, and these conditions are set to power the electrodynamic machine. 

Next, the voltage and current signals at the input of the electrodynamic machine 
are determined by direct measurement, by knowledge of imposed input conditions or 
other conventional means known in the art. 

Thereafter, positional and force parameters of the load device are estimated in a 
sensorless manner, from the system model and the current and voltage signals. It is 
appreciated that such estimating may be facilitated with system modeling via 
experimentation, mathematical methods, and/or combinations thereof. The system model 
only requires knowledge of the electrical and mechanical parameters of the linear 
machine , and not the load to which it transmits power. 
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The positional parameter to be estimated through sensorless means is preferably 
selected from the group consisting of displacement, velocity, and acceleration; and the 
force parameter is selected from the group consisting of force and pressure. As a group, 
these estimated parameters will be referred to as dynamical parameters. 
5 The estimated dynamical parameters are used either alone or in conceit with the 

obtained current and voltage and/or frequency, to determine at least one operating 
condition of the thermoacoustic device. The operating condition determined is an 
instantaneous operating condition that varies over time and with respect to changes in the 
load device. 

10 Next, the process continues with determining a difference between the estimated 

operating condition and a desired operating condition. It is appreciated that the desired 
operating condition for the thermoacoustic device may be one selected from the group 
consisting of maximum efficiency, maximum power and phase angle between the 
estimated dynamical parameters, as well as other operating conditions. In the case of 

15 maximum efficiency, or maximum power, the desired operating condition may not be 
known a priori, but may be identified via a searching through the input conditions. 
Finally, the voltage, current or frequency input to the electrodynamic machine is adjusted 
so as to reduce the difference between the instantaneous operating condition and the 
desired operating condition. In the case of maximum power or efficiency, the desired 

20 operating condition is the maximum value, so the inputs are adjusted to seek the 
maximum value, generally within operating limits of the device. 

BRIEF DESCRIPTION OF THE DRAWINGS 
To further describe the nature and objects of the present invention, reference 
should be had to the following detailed description taken in conjunction with the 
25 accompanying drawings in which like parts are given like reference numerals and 
wherein: 
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Figure 1 is a cross-sectional view of an exemplary thermoacoustic device which 
may be used with the present invention; 

Figure 2 is a cross-sectional view of another exemplary thermoacoustic device 
which may be used with the present invention; 
5 Figure 3 is a sectioned perspective view of the thermoacoustic device of Figure 2; 

Figure 4 illustrates an electrical equivalent model of a thermoacoustic system; 

Figure 5 generally illustrates a block diagram of a tracking circuit for tracking the 
electrical acoustical frequency of a thermoacoustic device wherein the tracking circuit 
comprises a phase lock loop controller; 
10 Figure 6 illustrates a comparison chart of estimated acoustic pressure and sensor 

obtained values for the same; 

Figure 7 illustrates a comparison chart of the estimated displacement and sensor 
obtained values for the same; 

Figure 8 illustrates a process flow block diagram for sensorless control of a 
15 thermoacoustic device as according to the invention; 

Figure 9 illustrates a flow chart of one embodiment of a method of sensorless 
control according to the present invention; and 

Figure 10 illustrates a flow chart of another embodiment specifically focused on 
maximizing efficiency or power. 

20 DETAILED DESCRIPTION OF THE INVENTION 

The present invention provides a method of sensorless control of a linear 
reciprocating electrodynamic machine used for driving a thermoacoustic device that 
includes a linear electrodynamic machine and a load device, such as a thermoacoustic 
refrigerator. However, it is appreciated that the sensorless control method described 
25 herein may be used with other devices with frequency dependent loads. 

Generally, the performance of the thermoacoustic device is controlled by 
obtaining an estimation of at least one dynamical parameter which has a mathematical 
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relationship to a desirable operating condition for the thermoacoustic device. Thereafter, 
the operating condition of the thermoacoustic device can be optimized by manipulating at 
least one input condition to the linear electrodynamic machine to achieve a desired 
operating condition relative to the estimated dynamical parameter. 
5 The invention provides a sensorless control approach for a thermoacoustic device 

where the acoustic pressure and actuator displacement at the load of the system are 
accurately estimated by using a pre-developed system model and the voltage and current 
at the input terminals of the linear machine. 

The exemplary embodiment, which we will call a first type of thermoacoustic 

10 device, herein relies on an operating condition at the thermoacoustic device in which it is 
desirable to maintain a relationship such that the acoustic pressure and actuator 
displacement are 90 degrees out of phase, to maintain an acoustic resonance. In another 
type of thermoacoustic device, which we may refer to as a second type of thermoacoustic 
device, a pure acoustic resonance is not present at the operating frequencies of interest, 

15 and it may rather be desirable to identify the condition of maximum efficiency of power 
delivery to the load. This condition is satisfied, when the ratio of the in-phase product of 
force and velocity (i.e. sensorlessly obtained dynamical parameters), divided by the in- 
phase product of current and voltage, is maximized, by changing the input conditions to 
the motor. Other types of thermoacoustic devices exist, which may have these or other 

20 desirable operating conditions as well. 

Typically for thermoacoustic devices, measurements from pressure sensors, 
displacement sensors, accelerometers or other sensing devices are used to determine the 
dynamic operating states of system parameters. Thereafter, control methods are used to 
adjust one or more of the input parameters to the linear motor to achieve a desired 

25 operating point at the thermoacoustic device. Hence, it is appreciated that one can effect 
the sensorless control of a thermoacoustic device to obtain a desirable operating point by 
providing a method for accurately estimating these dynamical operating parameters and 
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thereafter using the estimated parameters to make adjustments at the input of the 
thermoacoustic device to obtain a desired operating point. 

A variety of thermoacoustic devices are known in the art. An exemplary 
embodiment of a thermoacoustic device that can be used with the present invention is 
shown in Figure 1. This device is an example of the first type of thermoacoustic device 
mentioned above. Figure 1 is taken from U.S. Patent No. 5,647,216, the entire contents 
of which are incorporated herein by reference. A better understanding of this device may 
be had by reference to the incorporated patent. As shown, this thermoacoustic device is a 
double-ended device with drivers located at each end. Thermal components are located 
adjacent each driver, with the thermal components, including a stack, and a pair of heat 
exchangers at each end. The present invention may be used to control a device such as 
this, as well as other thermoacoustic devices. 

An example of the second type of thermoacoustic device discussed above is 
shown in Figure 2 and 3. Figure 2 shows a cross-sectional side view of the 
thermoacoustic device 110. Figure 3 shows the same thermoacoustic device 110 with 
portions cut away, and some components left out of the view for simplicity. The 
thermoacoustic device of Figures 2 and 3 is the subject of pending US Patent Application 
Serial No. 10/410,492, the entire contents of which are incorporated herein by reference. 
A more complete understanding of this device may be had by reference to the 
incorporated patent application. Additional thermoacoustic devices are discussed in 
pending US Patent Application Serial No. 10/409,855, the entire contents of which are 
incorporated herein by reference. The present invention may also be used to control one 
of these devices. 

The thermoacoustic device shown in Figures 2 and 3 will now be discussed in 
more detail, to give an understanding of this type of thermoacoustic device. The device 
110 is designed to produce a refrigeration effect by absorbing heat through a unitized 
cold head heat exchanger 112 and delivering that heat at a lower temperature to a 
regenerator 114, contained in a thermally-insulated and thermally insulating support or 
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platform 118. The energy in sound waves is used to pump that heat up to a higher 
temperature at the hot end of the regenerator 114. This heat, plus any acoustic energy 
deposited by the thermoacoustic heat pumping process, is deposited on the hot heat 
exchanger 1 16. The heat deposited on the hot heat exchanger 1 16 is exhausted from the 
5 system by hot heat exchange transport fluid that passes through the heat exchanger 1 16. 

The components of the thermoacoustic device 110 are preferably housed in a 
pressure vessel 120. The pressure vessel 120 is formed by a generally cylindrical wall or 
shell 122 that is closed on its "lower" end by a base plate 124. It should be noted that 
descriptors such as "upper" and "lower" are used merely for ease of description of the 

10 figures, but are not limiting on the configuration or positioning of this device. In fact, the 
device Figures 2 and 3 is preferably operated in an inverted position as compared to the 
illustrated position. It can also be oriented in other positions. The position descriptors 
used herein refer to the illustrated orientations. The upper end of the wall or shell 122 is 
closed off by the combination of the platform or support 118 and the cold head heat 

15 exchanger 112, as will become more clear with respect to the further description of the 
cold head heat exchanger 112 herein below. 

A linear motor 126 is attached to the bottom plate 124 and the moving portion of 
the linear motor 126 is connected to the main piston or power piston 128 by a connecting 
member 130. A bellows 132 extends between the power piston 128 and the support 118. 

20 This defines part of a compliant enclosure, in accordance with the disclosure in U.S. 
provisional patent application Serial No. 60/372,008, the entire contents of which are 
incorporated herein by reference. As the power piston 128 moves upwardly, the working 
fluid in the compliant enclosure is compressed, and as the power piston 128 moves 
downwardly, the working fluid in the compliant enclosure is decompressed. 

25 A multiplier volume 134 is defined inside the compliant enclosure. This 

multiplier volume 134 is defined by a multiplier cylinder 136, which has an upper end 
attached to the support 118 and a lower end closed off by a multiplier cone or piston 138. 
A flexible seal 140 interconnects the cone or piston 138 with the cylinder 136. As the 
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multiplier cone 138 moves upwardly and downwardly, the volume of the multiplier 
volume, and the pressure of the gas contained therein, is increased and decreased. 

The volume of working fluid inside the compliant enclosure and outside the 
multiplier volume 134 may be referred to as the main volume 144. The main volume 144 
and the multiplier volume 134 are in fluid communication through the thermal 
components. That is, a fluid path is defined from the multiplier volume 134, through the 
hot heat exchanger 1 16, through the regenerator 1 14, through the fins 142 of the cold heat 
exchanger 112, and through windows 146 in the support 118. Therefore, gas can flow 
from the main volume 144 through the windows 146, make a U-turn through the fins 142, 
and enter the regenerator 114. In operation, the multiplier piston 138 and power piston 
128 move generally in phase or close to being in phase with one another. Therefore, flow 
of gas between the multiplier volume 134 and main volume 144 is limited, since a 
pressure wave converges and diverges from both "sides" of the regenerator 114. As will 
clear to those of skill in the art, the device 110 allows for efficient refrigeration. 

Turning again to the present invention, the first step in performing the inventive 
method is to capture the pertinent electrical, mechanical, and acoustical dynamics of the 
thermoacoustic device for developing a system model. The thermoacoustic device may 
be any of those just described, or other devices now known or yet to be developed. In 
accomplishing this step, particular terms may be included or neglected in the model that 
are found to be unique to the linear machine comprised in the thermoacoustic device. 
However, it is appreciated that the methodology used to obtain the system model is 
intended to remain general. 

Electrical Model 
The electrical dynamics of the linear machine can be expressed as: 
.„ dA(xJ) 

where v, is the voltage excitation at the terminals of the machine, i is the stator current, 
R e is the stator winding resistance, and X is the flux linkage. Notice in the above equation 
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that core losses such as hysteresis and eddy current losses have been omitted. For a given 
core material, hysteresis losses typically vary linearly with frequency while eddy losses 
usually vary with frequency squared. Also, the winding resistance R e can reflect 
frequency dependent losses due to skin and proximity effects or other effects and can be 
determined empirically. Frequency dependent resistance losses can be added to the 
model if their effect is determined to be significant. 

The flux linkage of the linear machine could be a function of displacement, 
current, and even air gap. The radial air gap for the linear-electrodynamic machine 
considered here is relatively constant and it is not expected to impact the flux linkage and 
can be neglected. Assuming linear magnetic materials, saturation of the core is 
neglected, but this too can be modeled or evaluated empirically when deemed significant. 

The total flux linkage can be expressed as the flux linkage due to the self 
inductance of the windings L e and the permanent magnet flux-linkage Ap m : 

MxJ) = L e (x)i + A pm (x). 
The machine's back emf dA(xJ)/dt is described below as: 



dt e dt { dx dx 



dt 



The position-dependent component of magnetic co-energy W fld can be expressed as: 

= \L e (x)idi+ \\ m (x)di 
= U e (x)i 2 + A(x)i. 



Hence, the interaction force f/i d , between the permanent magnet field and the stator 
magnetomotive force (mmf) is: 
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20 



The first term mfju is the reluctance force. In a properly designed machine the reluctance 
force should be small, as the inductance is assumed to be independent of actuator 
displacement. 



dx ~ ' 



With reluctance forces neglected, fa becomes, 



. Mm**) ■ 
^ dx~ l - 



The change in permanent magnet flux linkage with respect to displacement is often 
assumed to be constant and referred to as Bl. 



Bl = 



v dx J 



10 The expression for the magnetic force applied to the actuator therefore simply becomes: 

Using /Lp m , dA(xJ)/dt, Bl and it follows that the electrical dynamics in 

.„ dA(x,i) 

v, = iR e + — , can be expressed as: 

at 

•n t di dx 
v t =iR+L— + Bl — . 
e e dt dt 

15 Mechanical Model 

The mechanical dynamics of the driver, bellows, and spring configuration can be 
modeled accurately with the traditional single degree of freedom mass, spring, and 
damper system. The spring constant is represented by K m or, alternatively, compliance 
Cm -1/K m , actuator moving mass as M m and the mechanical damping is represented by a 
dashpot with damping constant R m . Core losses, due to the interaction of the permanent 
magnet flux with the stator iron, can be modeled as part of the mechanical damping. The 
mechanical and acoustical dynamics are coupled by an actuating piston having effective 
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area A. The result is an acoustical pressure p a in the resonator at the piston location, with 
an attendant reaction force on the piston equal to the product paA- The velocity 
v m =dx/dt, with which the piston moves corresponds to a volume velocity 
u a =Av m -Adxldt in the acoustical chamber. The forces related to the spring, mass, 
damper, and external acoustic action are: 



fm ~M m -, 

J m m ^2 ' 



fa=Ap a , 

where x 0 is the equilibrium value of x with the spring unstretched. The spring is assumed 
to be linear and is operated around the equilibrium point of xo = 0. The mechanical 
dynamics can be represented as the sum of the forces: 

f M =Bli = M m ^ + R m ^- + K m x + Ap a . 

The acoustical impedance is -pju a , and when referred to the mechanical mobility 

domain using the effective piston area of the bellows, the acoustical system can be 
represented as: 

_1 

Z acs A 2 * 

Figure 4 shows an electrical equivalent model that represents the thermoacoustic system. 

Estimating State Variables 
The dynamic states of acoustic pressure and actuator displacement (the dynamical 
parameters) are needed, in the case of the first type of thermoacoustic device, to track the 
acoustical resonance with a PLL (phase locked loop). However, as noted earlier, it is 
appreciated that one may desire to effect a different optimum operating condition for the 
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system which would require different and/or other parameters (such as input conditions at 
the linear motor), in addition to a different control means. An example of the just 
mentioned kind, as noted earlier, applies to a thermoacoustic device of the second type, 
wherein one may desire to control for maximum efficiency, and this requires the 
dynamical parameters force and velocity, and the motor input voltage and current. 
Therefore, in general, the inventive method requires that at least one dynamical parameter 
be selected from the model for use in determining the instantaneous operating condition 
of the system. Thus, in the exemplary embodiment, acoustic pressure and actuator 
displacement are selected as the dynamical parameters. The next step in the inventive 
method for sensorless control of the linear machine involves accomplishing accurate and 
reliable estimations of these dynamical parameters. 

The selected dynamical parameters can be estimated using the physical model and 
knowledge of the input terminal voltage and current. In this case, the input terminal 
voltage and current signals represent parameters at the linear machine that have a 
functional relationship to the pressure and displacement, which depends on the load 
condition, that allows for accurate and reliable estimations to be obtained for the latter. In 
particular, in the case of maintaining an acoustic resonance, a known phase relationship 
between the pressure and displacement is desired, thus the input frequency is can be 
adjusted to achieve this. 

Thus per the general inventive method, the at least one input condition is set to an initial 
value and thereafter is adjusted to obtain the desired optimum operating condition in view 
of the estimated dynamical parameters. 

Hence, by manipulating the equation v =iR +L — + BI— wherein the time 

dt dt 

varying voltage and current represent the at least one input condition, the displacement 
can be estimated as: 




Once displacement is known, using the equation: 
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f JU =Bli = M m ^- + R m ^ + K m x + Ap a ., 
pressure can be estimated as: 

The electro-mechanical parameters R e , L et Bl y M m Rm, and K m can be found 
experimentally with least-squared curve fitting utilizing measured values of v,, z, and p a 
over the expected operating frequency range, and in general can be functions of 
amplitude and frequency. 

Actual signals are non-ideal and have small amounts of high frequency noise, and 
inherent dc offsets contribute to various circuit errors. The latter is especially true when 
the motor excitation is only approximately sinusoidal, as when driven by a pulse-width 
modulated signal. This can make the estimated pressure and displacement difficult to 
implement directly. If a signal with a dc component is integrated, the output of the 
integrator is unbounded, and likewise differentiators yield similar results for input signals 
with high frequency components. Since the frequency range of operation of a device 
controlled using the present invention is typically small, band-pass filters are used to 
attenuate all frequencies outside the narrow expected operating range. The filters delay 
the signal, which results in a phase shift, but care is taken to use the same filter for 
current and voltage reference signals during comparisons, thus producing an equal delay 
in all signals. If one signal does not require filtering, it is delayed to match the delay 
imposed by filtering the other signal. 

Voltage Frequency Control 

As discussed previously, it is desirable that the device described in this exemplary 
embodiment be operated at the acoustical resonance frequency. 

Typically for thermoacoustic devices, measurements from pressure sensors, 
displacement sensors, and accelerometers are used to determine the time varying state of 
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these parameters. Based on this parametric information, the voltage excitation v des to 
the linear machine can be adjusted until the desired states are achieved. To excite the 
acoustical resonance of the system, the frequency of the excitation voltage v des at the 
linear machine must be the same as the acoustical resonance frequency. It is appreciated 
5 that during normal operation of a thermoacoustic refrigerator the acoustical resonance 
frequency can change due to temperature changes in the thermodynamic working fluid 
gas, changes in the cooling load, or with changes in the inlet and outlet temperatures to 
cause the acoustic resonance frequency to depart from the excitation voltage frequency. 
As a result, the excitation voltage must be adjusted and made equal to the acoustic 

10 resonance frequency. Thus, the next step in the inventive method involves providing a 
control device capable of recognizing a difference between the instantaneous operating 
condition and a desired operating condition, and adjusting the excitation voltage 
frequency at the linear machine such that the desired operating condition (acoustical 
resonance frequency) of the linear machine is accurately and efficiently obtained. 

15 At the acoustical resonance, the pressure at the piston face and the volume 

velocity of the gas are in-phase. Hence, the pressure and displacement are 90 degrees out 
of phase. Since it is desired that pressure and displacement be 90 degrees out of phase, 
the simplest and arguably the most effective controller is the PLL (Phase-Locked-Loop) 
controller. This controller only provides the frequency, as the amplitude of the VCO 

20 (Voltage-Controlled Oscillator) output is controlled separately. The magnitude of the 
desired voltage \\v des \\ is dependent on the cooling load temperature and not the 

operating frequency, and represent a second input condition which could also be varied, 
depending on the desired operating condition. 

Figure 5 provides a schematic representation of usage of the estimated 
25 displacement x and the estimated pressure p a to drive the frequency of the excitation 
voltage to the acoustical resonance frequency. In Figure 5, LPF represents a low pass 
filter, which is used to select the time averaged component from the product of the 
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estimated displacement and pressure, and the "PI Controller" represents a standard 
proportional-integral controller, used to provide a control signal to the voltage controlled 
oscillator (VCO). The VCO controls the system excitation frequency, thus closing the 
control loop since the latter impacts the relationship between x and p a . 
5 Preferably, the method of sensorless control as according to the invention may be 

used to obtain a plurality of operating conditions in view of the dynamical parameters. 
Such conditions may illustratively include maximum power at the load, maximum 
pressure or force, or a desired phase relationship between dynamical parameters. 

Illustration of Performance 

10 As explained in the 'detailed description of the invention', to improve the 

agreement between measured and estimated values, the implementation of the position 
estimation can be modified based on the empirical observations, to include higher order 
terms in current that may be associated with non-linearity in the magnetic materials, and 
the frequency dependent resistance mentioned earlier. Such an implementation was used 

15 on a linear motor very similar to that described in "Linear Electrodynamic Machine" US 
Patent No. 5,389,844, to operate a thermoacoustic device (of the first type) described in 
the exemplary embodiment, wherein tracking of the acoustic resonance was desired. So, 
the modified position estimation for Figure 7 below is: 

* = i/[ Jl v ' " iR < " " w) * - V + *i' 3 + */'] . 

20 with ki and k 2 representing small, non-complex constants. The second modification 
observable in comparison with the estimated position is the explicit inclusion of a linear 
frequency dependence in the effective coil resistance. In the process of parameter 
extraction, the coefficient of the linear frequency dependent resistance appears as the 
imaginary component of the inductance, so the additional frequency dependant resistance 

25 is here represented by the product of L^g, and co, the operating frequency in radians/sec. 

The graphs in Figure 6 provide a comparison of the estimation of pressure at the 

face of the transducer, with simultaneously measured values, taken during the process of 
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parameter evaluation. Similar graphs are given for displacement in Figure 7. This level 
of agreement was found suitable for the operation of a thermoacoustic refrigerator of the 
type in the exemplary embodiment, but it is appreciated that other terms could be 
included, were this to be found necessary for a particular implementation of the 
5 thermoacoustic refrigerator, and selection of type of operating condition. 

Figure 8 provides an overview of a process flow diagram for the method of 
sensorless control of the thermoacoustic device as according to the invention. 

Figure 9 generally summarizes a method for sensorless control of a 
thermoacoustic device, and depicts the method as a flowchart. As shown, the first step is 

10 to set the electrical inputs to the electrodynamic machine. This typically means setting 
the voltage or current, and the frequency. Typically voltage and frequency are set, and 
current responds in a manner governed by the characteristics of the load and the 
parameters of the electrodynamic machine. Next, a current and voltage signal are 
determined at the inputs of the electrodynamic machine. One of these signals usually 

15 may be determined directly from the power supply, since it was set, while the other signal 
is measured. Alternatively both may be measured. The voltage and current signal are 
used to estimate a positional parameter, such as displacement. This is typically done 
using a mathematical model of the machine. The current signal and the positional 
parameter are then used to estimate a force parameter, such as pressure. An operating 

20 condition is then determined using at least the positional and force parameters. In one 
example, the operating condition is the phase relationship between displacement and 
pressure. This phase angle is proportional to the cycle averaged product of these two 
parameters. If they are 90 degrees out of phase, the product will be zero. In this same 
example, the desired operating condition is that the displacement and pressure are 90 

25 degrees out of phase, which corresponds to a zero product. A difference between the 
desired operating condition, a zero product, and the actual operating condition, is then 
determined, and the frequency, current or frequency are adjusted to decrease the 
difference. The process is repeated as necessary. For some situations, the operating 
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condition may depend on additional parameters, such as the electrical inputs, as well as 
the force and positional parameters. 

Figure 10 provides a flow chart for a method more specific to maximizing 
efficiency or power. This is similar to the above example, except that operating condition 
is power or efficiency and the inputs are adjusted to maximize the power or efficiency, 
typically within certain performance limits (maximum current, voltage, displacement, 
etc.). Determination of the power may be made based on the force and positional 
parameter. It may be calculated by multiplying force times velocity and time averaging 
this product over a cycle. Knowledge of the electrical inputs may be desirable to avoid 
exceeding the current or voltage capabilities of the machine. Also, maximum 
displacement, or other positional or force parameters may have to be limited to avoid 
damage. Determination of efficiency requires using the electrical inputs in addition to the 
force and positional parameters. The efficiency is the power to the load (product of force 
and velocity) divided by the electrical input power (product of voltage and current time 
averaged over a cycle). Again, limits may be placed on force or positional parameters, or 
electrical inputs. 

For cases wherein it is desirable to calculate maximum power or maximum 
efficiency, and for cases in which the load is something other than a thermoacoustic 
device, optimum selection of the dynamical parameters to be estimated, may be different 
than those described explicitly in the equations above. For purposes of illustration, it 
may be desirable to select among the dynamical parameters to estimate force and 
velocity, for operating points which are composed of these parameters. In such a case an 
example of an expression to estimate the velocity, V pist , is given by: 




and an example of the force on the load, f load , is given by: 
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f load =Bli-M m -^-R m V pisl -K$V pist dt 
where the parameters, are as noted above. 

The control techniques described herein can be applied to permanent magnet, 
moving coil, and other reciprocating linear electrodynamic drive machines and has utility 
5 for cases in which it is desirable to estimate the magnitude, and phase relationship 
between the force and displacement, as well as to adjust the relationship by means of 
adjusting the operating frequency. It is appreciated that this approach is not limited to 
thermoacoustic refrigeration. 

It should also be noted, that while the description herein refers to an input to the 

10 linear machine which has a single frequency (i.e. a pure sinusoid) the method is 
sufficiently general, such that the input excitation to the linear motor may be only 
approximately sinusoidal (i.e. a pulse width modulated signal, typical of the 
commercially available systems for controlling input frequency to rotary motors) as long 
as the current and voltage are suitably filtered or conditioned, for use in estimating the 

1 5 dynamical parameters . 

Another example of an application in which it is desirable to track an acoustic 
resonance is in the field of sonic gas analysis, wherein one can make use of the fact that 
the sound speed in a known gas changes with the addition of contaminants, and that when 
the gas is contained within a fixed wavelength resonator, the sound speed can be obtained 

20 from the resonance frequency. In the past this tracking has been accomplished with 
sensors, or separate microphones within the resonator, which this invention eliminates. 
In this case, the directly tracked frequency can provide information regarding changes in 
the mean molecular weight of constituents within the acoustic resonator. The present 
invention may more broadly be considered as a method for sensorlessly controlling a 

25 system including a reciprocating linear electrodynamic machine that harmonically drives 
a displaceable portion of a load device. With the exemplary thermoacoustic device, the 
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electrodynamic machine is a linear motor, the load device is the thermoacoustic 
refrigeration portion, and the displaceable portion is the piston. 

It is also possible to implement the position estimation in a phasor form. In such 
case, two harmonic signals at the commanded frequency are generated, with a 90-degree 
5 phase relationship between the two. Implementation of the equations then would be 
possible without use of numerical or analog integration and differentiation, and explicit 
reference to complex quantities. 

It is also clear that one could obtain the piston position via means other than the 
estimation used here and still apply the methods herein to obtain the magnitude of the 

10 force or pressure only, in order to provide information suitable for control of the linear 
electrodynamic machine. In this situation, the method as shown in Figure 9 would be 
modified such that the "estimate a positional parameter" step is replaced with a "measure 
a positional parameter step". The step of determining the voltage signal may become 
unnecessary since the force parameter may be estimated based on the measured, or 

15 otherwise determined, positional parameter and current. The voltage signal may be used 
for other reasons, as needed. This version of the present invention would still provide 
benefits since the force parameter could be estimated instead of measured, thereby 
eliminating one sensor. 

It is recognized that for particular control applications, it may be desirable to use 

20 operating points other than the phase between the force and positional parameters, or 
efficiency, or maximum power. In some applications it may be desirable, for example, to 
have the ratio of a force parameter and a positional parameter be defined as the operating 
condition. As a specific case of the latter, where the ratio is force divided by velocity, 
this operating condition would be the mechanical impedance. Control could be 

25 implemented, to tune the electrical excitation, so as to achieve a particular value of 
impedance, to reach a maximum value, or a minimum value. It is further recognized, that 
other operating conditions can be mathematically defined, which may be informed by the 
value of a force parameter and a positional parameter, and sensorless estimation of these 
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parameters permits control to any such operating point, assuming it lies within physical 
possibilities of the combined system of the machine and the load. 

While the invention has been particularly shown and described with reference to 
preferred embodiments thereof, it will be understood by those skilled in the art that 
5 various alterations in form and detail may be made therein without departing from the 
spirit and scope of the invention. 

We claim: 
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